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Statistical characterisation of full-disk EUV/XUV solar 
irradiance and correlation with solar activity 
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Abstract. We investigate the distribution of fluctuations in solar irradiance when integrated over the full disk, 
obtained using extreme ultraviolet/soft X-ray observations from the SOHO CELIAS/SEM instrument. This time 
series sums over both the contributions of single distinguishable flares, and of many other processes. By detrending 
we select events with timescales of less than a few hours such as waves, slow flows, and CMEs. The statistics 
generated by this range of phenomena can be characterised by power-law-tailed distributions. We show that (i) 
during the high activity period 2000 Jan-June the tail exponent ar = 1.5 ±0.1; (ii) during the low activity period 
1996 Jan-June ar ~ 3.0 ± 0.2; and (iii) in general ar decreases with increasing activity. 
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1. Introduction 

The statistical characteristics of solar irradiance have 
direct implications for coronal heating and terrestrial 
climate change; these statistics can be addressed both 
through flare events and, as in the present paper, through 
global measures. Many authors have identified transient 
flaring events and quantified their frequency distributions. 
The earliest observations of soft X-ray bursts by Drake 
H1974|l . and of hard X-ray bursts by Datlowe et al. I|1974|l . 
Lin et al. I|1984|l . and Dennis 119851, found that the num- 
ber of bursts with a given peak photon flux follows an 
inverse power law with logarithmic slope qe ~ 1.8. These 
findings gave rise to the model of the active corona as 
the superposition of a large number of impulsive phenom- 
ena, from the largest flares (10'^^ ergs) down to nanoflares 
{10^'^ ergs) and possibly beyond. The contribution of 
nanoflares to coronal heating is reviewed by Berghmans 
120021 Kucera et al. ( 1997) present evidence for changes in 
the upper limit of flare energies with the sizes of active 
regions, while secular variations are generally not seen. 
However, Bai l|1993|l and Bromund et al. H1995|l report a 
154-day periodicity in hard X-ray burst distributions, and 
Wheatland ( 2000J finds a variation in the fiaring rate. 

In this paper we investigate a different statistic: the 
distribution of fluctuations in the irradiance from the full 
solar disk, using EUV/XUV observations by the SOHO 
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CELIAS/SEM instrument (described in Sect. This 
global data encompasses both flaring and other activity 
such as waves, slow flows, and CMEs. Importantly, this 
also includes flare activity that could not be identifled by 
event selection from images. By characterising the distri- 
butions with single power-law tails, we find that the ex- 
ponent thereof differs by a factor of two between periods 
of high and low solar activity. 

2. SOHO/SEM Data 

The Charge, Element, and Isotope Analysis System/Solar 
Extreme Ultraviolet Monitor (CELIAS/SEM) transmis- 
sion grating spectrometer, on board the Solar and 
Heliospheric Observatory (SOHO) spacecraft, has mea- 
sured with 15 s resolution the full disk absolute photon 
flux at 1 AU since the launch of SOHO in 1995 Dec. FuU 
details of the instrumentation and calibration are given by 
Hove et al. (|1995|) and Judge et al. (1998j) . and the data 
are available via the internet^. Previous EUV/XUV mea- 
surements suffered from instrument degradation, but the 
SEM data is highly stable and reliable with an estimated 
absolute uncertainty of ~ ±10% (1 a) (Judge et al. l2UU2fl . 
Owing to the statistical correlation of electron tempera- 
ture and thermal energy in flaring processes (Aschwanden 
I1999f) . narrowband detectors provide biased estimates of 
flare energy distributions (Aschwanden & Charbonneau 

^ http:/ /www. usc.edu/dept/space_science/sem_first. htm 
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Fig. 1. Time-series of raw full-disk solar flux (observed 
at f5 s intervals by SOHO/SEM); (upper) 1996 Jan- June 
(863450 measurements), (lower) 2000 Jan- June (1028342 
measurements) . 



I2002|l : to overcome this problem, we choose the broad- 
band (0.1-50 nm) channel 2 data. In Fig. we plot the 
full-disk solar flux averaged over 15 s intervals from 1996 
Jan-June (a period of low activity) and 2000 Jan- June 
(a period of high activity). The six or so quasiperiodic 
variations with the largest amplitude are caused by the 
movement of active regions across the disk with the sun's 
27-day rotation period, plus contributions from the evolu- 
tion of active regions over days to months (Frohlich & Pap 
I1989|l and from the 1 1-year solar cycle (Willson & Hudson 

Since we wish to study only the short-time fluctua- 
tions, it is necessary to remove the long-term variations 
from the time-series. To achieve this we use a technique 
known as detrending, and the results are shown in Fig.|5] 
Detrending is effected by subtracting from each data point 
the mean of the measurements over a certain time r cen- 
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Fig. 2. Time-series of raw full-disk solar flux (observed at 
15 s intervals by SOHO/SEM) detrended by subtracting 
from each measurement the mean of the 10000 measure- 
ments (spanning r = 40 hours) centred in time on that 
point; (upper) 1996 Jan- June, (lower) 2000 Jan- June. 



tred on that point, t is chosen to be much less than 
the timescales to be removed and much greater than the 
timescales of interest. Here we set r = 40 hours (10000 
15 s sampling intervals) so as to remove the dominant 27- 
day cycle while retaining variations on timescales of up to 
a few hours. 



3. Full-disk irradiance distributions 

In Fig. 13 we plot the number distributions N{x) of the 
detrended intensities x during periods of low and high ac- 
tivity (shown as time-series in Fig.[2Jl, using log-linear axes 
to magnify the small probabilities in the tails of the dis- 
tributions. It proved impossible to fit single distribution 
functions to the entire range of detrended data in this fig- 
ure. Gaussian distributions, which are parabolic on these 



J. Greenhough et al.: Full-disk EUV/XUV solar irradiance and solar activity 



X 



^ 4 

cn10 

c 

o 




2 4 6 

Intensity x (10®photons cm"^s"^) 




1 2 

Intensity x (10^°photons cm"^s~^) 



Fig. 3. Log-linear distributions logA^(a;) of detrended 
full-disk solar flux x (observed at 15 s intervals by 
SOHO/SEM); (upper) 1996 Jan- June, (lower) 2000 Jan- 
June. Vertical error bars indicate a factor of ±1/ y^N{x) 
for bins having N{x) > 2; bin widths (not shown for 
N{x) = 0) indicate horizontal errors. Dashed lines indi- 
cate the lower limits for the fits in Fig. 01 

axes, clearly do not apply; log-normal distributions are 
inappropriate since many detrended values are negative. 
Frechet distributions (Chapman et al. 120021 Greenhough 
et al. nna^ may be fitted but they fit only ~ 30% of the 
bins within the errors shown in Fig. 3. Figure 4, which 
plots on log-log axes the ranges to the right of the dashed 
lines in Fig. 3, shows that the tails of the distributions 
approximate well to power laws. Single power laws of the 
form N(x) ~ where ot is the tail exponent, can 

be fitted to give (i) or = 3.0 ± 0.2 for 1996 Jan-June, 
and (ii) ut = 1.5 ± 0.1 for 2000 Jan- June, both having 
an goodness-of-fit of 0.99. To obtain the correct value 
for qt from the slope m of logiV(a;) against log a: (since 
Ot — ~n^), the data are placed into bins whose widths 
increase exponentially with x. One may also fit composite 
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Fig. 4. Log-log tails of distributions log N{x) of detrended 
full-disk solar flux logx (observed at 15 s intervals by 
SOHO/SEM); (upper) 1996 Jan-June, (lower) 200 Jan- 
June. Vertical error bars indicate a factor of ±l/^/N{x); 
bin widths indicate horizontal errors. Slopes m are calcu- 
lated via least-squares fitting with 95% confidence inter- 
vals; tail exponent ax = —m. 



functions such as broken power laws. If the fitted slopes 
are allowed to change at, for example, x — 10^ (upper 
plot) and a; = 5 X 10^ (lower plot), the best fits are re- 
spectively ar = 2.8 ± 0.2 ^ 3.3 ± 1.2 (overall = 0.95) 
and OT = 1.3 ± 0.2 1.5 ± 0.2 (overall i?^ = 0.97). 
Importantly, regardless of the details of the chosen fits, we 
find a significant difference in the tails of the detrended 
full-disk solar irradiance observations between the selected 
six-month periods of low and high activity. 

In Fig. 13 we use single power laws as the simplest de- 
scription of the distribution tails for intermediate times. 
We plot the slope —ar (with error bars showing 95% least- 
squares confidence bands) for each available month of data 
from 1996 Jan to 2000 Dec, against the sunspot number 
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Fig. 5. Slopes m — —ar of log-log distributions of de- 
trended full-disk solar flux for each available month of data 
1996 Jan-2000 Dec, plotted against the sunspot number 
for that month. Vertical error bars indicate 95% confidence 
intervals for the least-squares fits used to obtain the val- 
ues of m; three values whose errors span a greater range 
than the plot are omitted. 

for that month (a proxy for solar activity)^. Each value is 
calculated using a least-squares fit from the bin proceed- 
ing the mode up to the last non-empty bin in the log-log 
distribution of detrended data for that month. In three 
cases - 1998 Oct, 1999 Feb, and 2000 Nov - the error bars 
extend beyond the range of all the calculated slopes, and 
these points are not included in the figure. The plot is 
suggestive of a trend, albeit with a large variance, toward 
flatter power-law slopes as the sunspot number increases, 
which is consistent with the six-month averages presented 
in Fig. 21 To our knowledge, these results represent the 
first clear activity-related distinction in the statistics of 
full-disk solar irradiation. 

4. Conclusions and discussion 

By removing the low-frequency (greater than a few hours) 
variations in full-disk EUV/XUV solar irradiance obser- 
vations, and calculating the distributions of the remaining 
detrended fluctuations, we obtain distributions whose tails 
can be characterised by single power laws. The power-law 
exponent shows significant variation with activity: for the 
low-activity period 1996 Jan- June, we find ot = 3.0 ±0.2, 
while qt = 1.5±0.1 for the high-activity period 2000 Jan- 
June. Comparing ax and the sunspot number for each 
month of data suggests that there is in general a flattening 
of the slope with increasing activity. These results provide 
evidence for activity-related changes in irradiance fluctu- 
ations of the full solar disk. This is particularly intriguing 

^ Available at|http://sidc.oma.be/lLtml/sunspot.htnil| 



given the diverse range of phenomena that contribute to 
this integrated emission measure. 

Changes in the value of ot may reflect important vari- 
ations in the proportion or character of emissions from 
non-flaring regions. However, uncertainty is introduced by 
the dependence of the absolute photon flux on the rela- 
tive spectral distribution, which may change during the 
intense flaring events (Ogawa et al. tHl98|i) that are more 
frequent during high activity. One would have to remove 
the XUV component from the EUV/XUV total - a pro- 
cess involving much uncertainty (Judge et al. I2002|l - in 
order to test for this effect . Furthermore, since the re- 
lationship between total quantities and their underlying 
components is an unsolved problem in the field of non- 
equilibrium statistical mechanics (discussed by Chapman 
et al. H2002|l . and by Vekstein & Jain H2003|l in the context 
of nanoflares), we do not attempt to make a quantitative 
link between a-r and the flare statistic ue (see Sect.^ in 
this paper. 
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